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Abstract 

We report two measurements of the top quark mass M top using the CDF II detector at the 
Fermilab Tevatron in a 318 pb -1 data sample observed in the lepton + jets final state. One method 
uses an event-based likelihood technique resulting in M top = 173.2 ^4 (stat.) ±3.2 (syst.) GeV/c 2 
or 173.2 GeV/c 2 . The second method reconstructs a top quark mass in each event using the 
measured invariant mass of the hadronically decaying W boson to constrain the jet energy scale 
to obtain a value for M top of 173.5 tH (stat.) ± 1.3 (syst.) GeV/c 2 or 173.5 ±g;| GeV/c 2 . We take 
the latter, which is more precise, as our result. 

PACS numbers: 14.65.Ha, 12.15.Ff 
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The top quark is the heaviest known elementary particle with a mass approximately 40 
times that of the next-heaviest quark or lepton. Because of this comparatively large mass, 
top quark studies provide insight into our understanding of mass in general, and test theories 
that explain the large range of quark and lepton masses. Within the context of the Standard 
Model of particle physics, the top quark mass is related to the masses of the W boson and the 
Higgs boson, the latter object being the key to our understanding of the origin of mass ]]. 
Precision measurements of the top quark and W boson masses test the consistency of the 
Standard Model, and in particular the Higgs mechanism. A precision measurement of the 
top quark mass is therefore a main goal of the experiments at the Fermilab Tevatron collider. 

In this Letter we present two measurements of the top quark mass in the lepton + jets 
decay channel. We use a sample of tt decays corresponding to 318 pb _1 of proton-antiproton 
collisions at \/s = 1.96 TeV and collected using the Collider Detector at Fermilab (CDF II) 
between February 2002 and August 2004. In the lepton + jets channel, tt pair production is 
followed by the decay of each top quark to a W boson and a b quark, the hadronic decay of one 
W boson, and the leptonic decay of the other. This decay channel has the largest branching 
fraction with good signal-to-background, allowing accurate top quark mass measurements. 
Events in this channel contain an electron or muon and a neutrino from the leptonic W 
boson decay, two quark jets from the hadronic W boson decay, and two 6-quark jets. 

We select events consistent with this decay topology and analyze them using two com- 
plementary methods. The first method uses an event-by-event likelihood analysis employing 
the leading order matrix element for tt production and decay to extract a joint likelihood 
as a function of the top quark mass, M top . This technique, known as the "dynamical likeli- 
hood method" or DLM, was developed by the CDF collaboration |2( and is similar to that 
used by the D0 collaboration to make the previous most precise measurement of the top 

n n 

quark mass |3j. The second method, developed by the CDF collaboration |4J], reconstructs 
a top quark mass, m£ eco , in each event and compares the distribution of m[ eco with tem- 
plate distributions derived from model calculations to estimate M top . We have improved 
this "template method" by making further use of the fact that the hadronically decaying W 
boson daughters should form a final state whose invariant mass is consistent with the known 
W boson mass and width. This allows us to constrain the jet energy scale, an important un- 
certainty in the earlier measurements. These two methods have a top quark mass accuracy 
30% greater than earlier results, and have different statistical and systematic uncertainties. 
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The CDF II detector |5J] is a general-purpose charged and neutral particle detector located 
at the Tevatron collider. We employ cylindrical coordinates where 9 and are the polar 
and azimuthal angles, respectively, with respect to the proton beam, and pseudorapidity is 
rj = — lntan(#/2). Transverse energy and momentum are Et = Esin8 and pt = psin.9, 
respectively, where E and p are energy and momentum. The detector comprises a solenoidal 
charged particle spectrometer, consisting of an eight-layer silicon microstrip detector array 
and a cylindrical drift chamber immersed in a 1.4 T magnetic field, a segmented sampling 
calorimeter with acceptance up to pseudorapidity \r)\ = 3.6, and a set of charged particle 
detectors outside the calorimeter used to identify muon candidates. 

Events for these analyses were selected by requiring an electron or muon candidate with 
transverse momentum p^ > 20 GeV/c and \r)\ < 1 and missing transverse energy exceeding 
20 GeV, corresponding to a high-energy neutrino candidate. The signal-to-background was 
improved by requiring in each event the presence of four or more jets with \r)\ < 2.0. To 
reduce backgrounds further we required either (a) at least four jets with transverse energy 
Et > 21 GeV or (b) at least three jets with Et > 15 GeV and a fourth jet with Et > 8 GeV 
with at least one jet with Et > 15 GeV identified as a b quark candidate through the presence 
of a displaced vertex within the jet arising from the decay of the long-lived bottom hadron 
(b tag). This selection resulted in 165 events that, based on our background estimates, are 
primarily ti events. The methods used to estimate the backgrounds are detailed in Q]. 

The DLM analysis uses a 63-event subset of those data defined by requiring exactly four 
jets with Et > 15 GeV where at least one of the jets has a b tag. We have estimated the 
various sources of background contamination in this sample, summarized in Table HI to be 
9.2 ± 1.8 events. The template method divides the 165 events into four non-overlapping 
subsamples with different expected m[ eco distributions and background levels. Ordered by 
decreasing statistical power, the subsamples are 1) events with at least four jets with E T > 
15 GeV and one 6-tagged jet ("1-tag Tight" sample with 63 events), 2) events with two 
or more 6-tagged jets ("2-tag" sample with 25 events), 3) events with a fourth jet with 
8 GeV < Et < 15 GeV and one 6-tagged jet with Et > 15 GeV ("1-tag Loose" sample with 
33 events), and 4) events with four jets with Et > 21 GeV and no 6-tagged jets ("0-tag" 
sample with 44 events). The estimated background levels in the samples with a b tag are 
summarized in Table UJ The background level in the 0-tag sample is determined in the 
subsequent fit. 
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TABLE I: The background composition and the number of ti candidates for events with > 1 b tag, 
and for the subset used in the DLM analysis. Statistical and systematic uncertainties have been 
combined. 





> 1 b tag 




DLM sample 


Source 




Expected Background 




W + jets 


19.6 ±2.4 




5.3 ± 1.1 


Multijet 


4.7 ±0.7 




3.1 ± 1.0 


Other 


2.3 ±0.2 




0.8 ±0.1 


Total 


26.6 ±3.0 




9.2 ±1.8 






Selected ti Candidates 




Data 


121 




63 



Both analyses use calibrated jet energies, based on a combination of instrumental cali- 
bration and analysis of data control samples []|. The uncertainty a c on the jet energy scale 
is the a major source of uncertainty on M top . For jets in the ti sample, a c is approximately 
3% of the measured jet energy, depending on the r] and pt of the jet. The parameter JES is 
defined as the difference, averaged over all jets in the sample, between the true jet energy 
and our measured jet energy after calibration in units of o c . 

The DLM technique, described in detail in j^] , defines a likelihood for each event based on 
the differential cross section per unit phase space volume of the final state partons, da t t/d&, 
as a function of M top . Detector resolution effects are accounted for using ti events generated 
by the HERWIG Monte Carlo program and full detector simulation to derive a transfer 
function (TF). The TF relates the transverse energies of the quarks, denoted by x, and the 
observed jets. For a given event, a Monte Carlo integration is performed over the possible ti 
final state kinematics in the following way: we first generate a random value for the virtual 
mass squared of the W boson in the leptonic channel, sw, according to the Breit-Wigner 
form. We identify the momentum of the electron or muon daughter with the measured 
value, and the neutrino transverse momentum with the measured missing transverse energy. 
We then generate random values for the momenta of final state quarks according to the TF 
probabilities. We determine the z component of the neutrino momentum, with a two-fold 
ambiguity, using sw as a constraint. Thus, for a given set of x and s\y, we fully determine 
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the event kinematics, and the event likelihood as a function of M top is given by 

L(M top ) = nJ2J2 ^jf^ x > Sw ^ w 

Ij Iu 

where the normalization factor N is independent of M top for a given event, and the indices 
Ij and I v run over the parton-jet assignments and the two neutrino solutions, respectively. 
The event likelihood is obtained by numerically integrating over x given by the TF and sw 
given by the Breit-Wigner distribution. 

Figure ^ shows the distribution of the top quark mass value at the point of maximum 
likelihood in each event compared with the expectation from simulated events. An inset 
shows the joint log-likelihood as a function of M top , formed by multiplying the likelihoods 
of the individual events together. We account for the presence of background events by 
evaluating the shift of —1.4 GeV/c 2 they make in the measured top quark mass. From the 
joint likelihood we infer M top = 173.2 1J; 4 (stat.) GeV/c 2 , where the uncertainty is only 
statistical. The systematic uncertainty due to the jet energy scale is estimated as the shift 
in M top arising from a 1 a c change in JES, and is 3.0 GeV/c 2 . 

The template method is described in detail in [lOj. We perform a y 2 minimization to fit 
the parton momenta from the tt daughters and determine m[ eco for each event, assuming 
that the final state arises from the decay of a tt pair into W bosons and b quarks. We use 
only the four leading jets in the mass reconstruction. In the x 2 fit> both sets of W decay 
daughters are constrained to have the invariant mass of the W boson, and both Wb states 
are constrained to have the same mass. The ambiguity arising from the different ways of 
assigning the jets to the four quarks is resolved by selecting the assignment with the lowest 
X 2 , taking into account the 6-tagging information. We construct a histogram of mj; eco for 
each subsample, discarding events with x 2 > 9, corresponding to poorly reconstructed or 
background events. 

The parameter JES is determined within this event sample by removing the W boson 
mass constraints, and identifying for each event all pairs of jets that would be consistent 
with the W boson final state. We form histograms of the invariant masses of these jet pairs 
for each of the four event subsamples and compare these with what we expect given the 
precisely known W boson mass 

We use these eight histograms to measure simultaneously M top and JES. An unbinned 
likelihood fit is performed to parameterized signal templates taken from simulated ti events 
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FIG. 1: The value of the top quark mass at the maximum of the DLM likelihood is plotted for 
each event. Data events (points) are compared to an expected distribution (histogram) comprising 
simulated ti (Mt op =172.5 GeV/c 2 ) and background events. The last bin includes events with 
masses > 305 GeV/c 2 . The inset shows the joint log- likelihood for the 63 events, before accounting 
for the presence of background. 

generated using different values of M top and JES, and background templates derived from 
studies of the relevant background processes. We include in the fit a Gaussian constraint 
(JES = Oil er c ) from the extrinsic jet energy calibrations, and we constrain the background 
rates in the 2-tag, 1-tag Tight, and 1-tag Loose samples to the estimated background rates 
within their uncertainties. The background level in the 0-tag sample is determined to be 
15.7 (stat.) events by the fit using the differences in predicted signal and background 
mass distributions. 

The four reconstructed top quark mass distributions and the results of the fit are 
shown in Fig. |21 where we also show the background contributions. In all cases, we 
see agreement between the observed data distributions and the fitted curves. We obtain 
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FIG. 2: The m^ eco distribution of the template method is shown for each subsample overlaid with 
the expected distribution using M top , JES and signal and background normalizations from the 
combined fit. The events with x 2 > 9 have been excluded from each subsample. 

M top = 173.5 jlg'g (stat.) GeV/c 2 , where the uncertainty is statistical and incorporates the 
uncertainty due to JES, which we estimate contributes ~ 2.5 GeV/c 2 . Figure El shows the 
likelihood in the M top -JES plane. If we do not constrain JES to the nominal value of zero, 
we obtain JES = —0.25 ± 1.22 cr c , which indicates our nominal jet energy calibrations are 
in good agreement with information provided by the W boson mass peak in the tt decay. 
This also demonstrates that the JES constraint from the W boson decay has comparable 
precision to the jet energy calibration. 

There are a number of additional systematic uncertainties that affect both analyses: 
initial state and final state radiation uncertainties (ISR/FSR), uncertainties arising from 
the parton distribution functions (PDFs), and uncertainties arising from modeling of the 
background processes, the choice of event generators and 6-jet fragmentation, decays, and 
color connections (Modeling) 12, Table ITU summarizes these uncertainties. 
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FIG. 3: Contours of the template method likelihood are shown in the M top -JES plane for the 
combined fit to all four subsamples. The crosshair shows the best fit point. Contours are given at 
intervals of AlnL, the change in log-likelihood from its maximum. 



TABLE II: The systematic uncertainties for the two analyses. 



Systematic 


DLM 


Template 




AM to p (GeV/c 2 ) 


AM top (GeV/c 2 ) 


Jet Energy Scale 


3.0 


[~ 2.5] a 


ISR/FSR 


0.6 


0.7 


PDFs 


0.5 


0.3 


Modeling 


0.7 


0.9 


Method 


0.5 


0.6 


Total 


3.2 


1.3 a 



a The JES uncertainty is included in the uncertainty reported by the likelihood fit. 



14 



The DLM method has additional uncertainties that arise from the use of transfer functions 
and from the procedure that corrects the measured mass for the presence of background 
(Method). Together with the JES and other common sources noted above, the systematic 
uncertainty on the DLM mass measurement is 3.2 GeV/c 2 . 

The template method has additional uncertainties arising from the statistical precision 
of the templates themselves and approximations made in treating JES as a single parameter 
affecting all jets coherently (Method). The total systematic uncertainty on the template 
mass measurement is 1.3 GeV/c 2 . 

In summary, we have presented two new measurements of the top quark mass. The anal- 
ysis using the DLM method results in M top = 173.2 + 2 ;^ (stat.) ± 3.2 (syst.) GeV/c 2 ; the 
analysis using the template technique results in M top = 173.5 (stat.)±1.3 (syst.) GeV/c 2 . 
There is a large statistical correlation between these measurements given the common 
data sample, so that we quote as a result only the more accurate measurement, the 
template method result of M top = 173.5 GeV/c 2 . This provides the most precise 
single measurement on this important physical parameter. In comparison, the previous 
most precise measurement was M top = 180.1 ± 5.3 GeV/c 2 [3] and the world average was 
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